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We present a continuum model for a nematic elastomer network formed in a chiral environment, 
for instance in the presence of a chiral solvent. When this environment is removed, the network 
can retain some memory of its chiral genesis. We predict the residual chiral order parameter for 
a number of possible scenarios, and go on to examine the robustness (stability) of the imprinted 
chirality. We show that a twist-untwist transition can take place, which determines whether the 
imprinting has been successful. A transition is via a coarsening of the helical director pattern and 
a lengthening of its pitch. Finally, the effect due to a subsequent swelling by an achiral solvent, or 
by a solvent of differing chirality, is considered. 



PACS numbers: 61.30.-v, 61.41. -fe, 78.20.Ek 



Nematic elastomers combine the properties of a liquid 
crystal and those of a conventional rubber. This synergy 
gives rise to novel material behaviour, which in turn has 
stimulated much research in past years y . From a gen- 
eral symmetry argument, de Gennes p| first suggested 
that chirality may be introduced to such an elastomer 
by simply forming it in a chiral solvent. The originally 
achiral liquid crystalline polymer would then remember 
the induced chirality after crosslinking, even when the 
solvent is removed or replaced with an achiral one. This 
is the case of chiral imprinting, which potentially can 
open up an entirely new way of producing materials of 
specified optical properties. Chiral imprinting, in prin- 
ciple, is akin to crosslinking a nematic polymer under 
an external magnetic or mechanical field where the 
monodomain state is also permanently imprinted. Ex- 
perimentally chiral imprinting was studied long ago 
as a function of solvent exchange. Recently imprinting 
has been studied as a function of both solvent re- 
moval and temperature. Another measure of imprinting 
occurs in intrinsically cholesteric networks. On temper- 
ature changes that would cause a substantial pitch vari- 
ation in a non-crosslinked cholesteric polymer melt, the 
corresponding network suffers essentially no variation - 
see for example Fig. 8 of reference where also many 
other references to cholesteric elastomers are given. In 
this paper, we analyse chiral imprinting, and predict the 
retained chiral properties of the elastomer when the ini- 
tial chiral environment of crosslinking is altered. Gradi- 
ents of director variation are modeled within continuum 
Frank nematic elasticity. The nematic elastomer penalty 
for rotation of the director relative to the solid matrix 
is described in a fully non-linear, rubber-elastic manner 
since rotation can be large. 

A nematic liquid crystal has a mobile director n, the 
gradient of which incurs a Frank energy 0| . For the free 
energy density of a cholesteric liquid crystal, the twist 
term is modified by a pitch wave number, q^: 



where i^i,2,3 respectively measure the energy penalty for 
splay, twist and bend, the three possible modes of the ne- 
matic director distortion. In a pure twisting, cholesteric 
conformation of the director, we can drop the terms in- 
volving Ki^j,. These other contributions arise if the di- 
rector moves away from the helical plane and tilts toward 
the pitch axis. The Frank energy can be viewed as a con- 
tinuum description with higher order spatial derivatives 
truncated, and it suffices for cases where director varies 
slowly over a nematic coherence length (~ lOnm). 

When liquid crystalline polymers are crosslinked into a 
rubber network, additional constraints on n arise in the 
form of director anchoring to the network. Anchoring 
manifests itself with an extra energy cost |^ for a uniform 
director rotation U3 relative to a local rotation of the elas- 
tic matrix fJ. The energy density is \Di [(Jl — u;) x n]^. 
A second term couples the relative director-matrix rota- 
tion to the shearing part of the elastic strain. A, that is 
D2 n . A . [{n — a;) X n]. For such strains to rotate the 
director in a cholesteric, they would have to vary along 
the helical axis. By the requirement of elastic compat- 
ibility this introduces secondary shears. One can show 
that these arc prohibitively expensive for this geometry. 
We accordingly ignore A terms. 




FIG. 1. A cholesteric director configuration. 

Consider an initially regular cholesteric helix along the 
z axis, Fig. 1. 6* is the azimuthal angle the local director 
makes within the x ~ y plane and is initially 
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where qo is the chiral pitch wavenumber. One can gener- 
alise the Di term to the large angle limit by considering 
a molecular nematic rubber elastic model It is no 
longer ^Di {9 — 6'o)^ but ^Di sir? {9 — q^z) which has 
the correct locally nematic symmetry of n = — n. The 
energy for an elastomer formed under a cholesteric sol- 
vent which is subsequently replaced with an achiral one 
is then: 

F = j dz l[K2 9'^ + Di sm\9 - qoz)] (2) 

where ' indicates d/dz. The first term in F wishes to 
remove the twist in n since now the chiral imperative 
of the solvent is removed the usual Frank twist penalty 
is fully incurred. However the second term insists that 
n is anchored to a helix thanks to the crosslinks being 
formed under cholesteric conditions. The initial network 
polymers are taken to be achiral, thus go is only induced 
and can be tuned with the choice of the chiral solvent we 
subject our elastomer to at formation. If achiral solvent 
is used when crosslinking, go could simply be zero, and 
we retrieve the description of more conventional nematic 
elastomers [0. 

The two limits of the energy density are (i) the per- 
fectly twisted cholesteric state / = ^K2qQ where the cur- 
rent pitch wavevector is unchanged from go, and (ii) the 
untwisted state / — jDi, the additional factor of ^ aris- 
ing from the averaging of sin^ over one period. Thus, 
crudely, we expect the director to be twisted if 

i^2go < Di/2 (3) 

and untwisted otherwise, this balance being tuneable 
since K2, Di and qg vary relatively to each other with 
temperature, degree of crosslinking and swelling and the 
presence of additional chiral agents. Eq. (||) anticipates 
the physics of our detailed results: highly twisted states 
{qo large) or systems with a large twist constant K2 
will pay a very high Frank penalty on loss of sponta- 
neous twist arising from the loss of the chiral solvent. A 
large combination K2qQ will overcome the anchoring Di 
and the elastomer will untwist - imprinting will be lost. 
Weakly twisted elastomers with weak twist constants will 
not overcome director anchoring and imprinting will re- 
main. We now analyze Eq. ^ for details of the phase 
behaviour. 

To simplify matters, we begin with the following sub- 
stitutions: 

(f) = qQZ-9 + Tr/2; u^z/^; £,^^/K^Jd^; a = ^go 

the angle 4> describes the variation away from (or modu- 
lation of) the original helical pattern, that is we are now 
in a rotating frame of reference. Lengths are reduced 
by the nematic rubber penetration depth ^, the natu- 
ral length scale in the problem. The parameter a is a 
non-dimensional measure of the nematic length relative 
to the chiral pitch, and the condition (pf) is equivalent 



to a '--^ 1/a/2- Following the remarks below Eq. (||), we 
expect imprinting to be lost at large a and retained at 
small a. Dropping a constant arising from the variable 
change 9 (f>, the reduced energy (per unit cross-section 
area perpendicular to the pitch axis) is: 

2F f 

- — = J du [{^~af~ sin^ (4) 

where (') signifies d/du. If we make the analogy to La- 
grangian dynamics, the integrand is the Lagrangian den- 
sity {L — T — V) oi a. particle in a potential given by 
sin^ 0. The problem also now resembles the problem of 
an electric or magnetic field applied perpendicularly to 
the helix of a liquid cholesteric, solved long ago by de 
Gennes and Meyer ||^. In reduced terms, the siv? (f) is 
like the field competing with the natural chirality a. In 
terms of the original problem, Eq. (||), it is as if a nat- 
urally untwisted nematic has a spatially chiral electric 
field "Z?i" applied to it in an attempt to induce a twist. 

The first integral of the Euler-Lagrange equation corre- 
sponding to Eq. (^ , the elliptic equation or often called 
in the literature the Sine-Gordon equation, leads to: 

0^ + sin^ = (5) 

where is the integration constant, which has the physi- 
cal interpretation of total energy. The analogy is helpful, 
and before we examine the details we already foresee two 
scenarios: If < 1, the particle does not have sufficient 
energy to overcome the barriers in sin^ (/) potential and 
is therefore localized, i.e. — 7r/2 < cf) < Tr/2. This corre- 
sponds to the case where cholesteric pitch is largely main- 
tained, with <^(m) only introducing small modulations to 
the director orientation. If > 1, our particle can climb 
out of the potential valleys and travel freely; this corre- 
sponds to the case of winding/unwinding the cholesteric 
twists in director orientation. Below, we determine c in 
terms of a. 

The localised limit, < 1. Our particle oscillates be- 
tween two values of (f)m — i arcsin c. We accordingly in- 
troduce a new variable (3 in the interval f3 G [—^, so 

that sini^ = c sin/3. Rewriting the derivative (j) in terms 
of $ and returning it to Eq. (0) reduces this equation to 

one of the standard elliptic form: /3 = \/l — sin^ p. 
The period of the oscillatory motion is found to be: 

Ti = 2 ^-2/C(c). (6) 

Here /C(c), and later £{c), are the complete elliptic inte- 
grals of the first and second kind respectively. The period 
Ti gives, in units of the characteristic length scale ^, the 
spatial repeat distance of our (/)-modulation of the origi- 
nal cholesteric angle qgz. The reduced energy of a period 
can be obtained from Eq. (^) as: 

Fi(Ti) = 2[(c2 + a^-2) /C(c) +2 £{c)]. (7) 
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In order to compared the stability of various states, we re- 
quire the reduced energy density. In the locahzed regime, 
< 1, we denote this by gi{c): 



g,=F,iTi)/Ti 



2 2 „ 2f(c) 



/C(c) 



(8) 



We now require the energy density in the traveling 
regime. 

The traveling limit > 1. The modulation period 
down the original helix corresponds to the time taken for 
our particle to travel from one peak to the next in the 
potential. It is calculated using Eq. (||): 



-r- = 2fc/C(fc) 



(9) 



where k = 1/c, and thus k < 1. The reduced energy of a 
period is: 



F2{T2) = 2[{a^ - c^)k IC{k) + 2c £{k) - a 
The corresponding energy density, (72(c) for > 1, is: 



(10) 



kJC{k) 



an 



(11) 



Combining (?i(c) and 52(c), Eq.s (j^, pT| ), we now have 
the energy density g{c) for the entire range of c. Fig. 1 
illustrates the energy density plots for three different val- 
ues of of the chiral strength, described by the parameter 
a. 
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FIG. 2. Reduced energy density, (7(c), for a = 0.2, 2/7r, 1. 
One has 51(c) for (? <\ and 52(c) for c^ > 1. 



We can see that at the transition point c = 1 the two 
energy densities approach each other. However, the den- 
sities are not smoothly joined, forming a cusp which turns 
from pointing downward to upward upon the parameter 
a increasing past a critical value ac = 2/7r. Invoking the 
identities for differentials of complete elliptic integrals: 

die £ IC d£ £ ~ K- 

Ik " fc(l-fc2) " 1 ' dk^ k 



we can minimise the reduced free energy with respect to 
k, or equivalently c, by setting dg/dc = 0. As evident 
from Fig. ^ minima only exist in the region c'^ > 1, that 
is in 52(c), and only when a > ac- There are no minima 
in gi. The analogous problem of a cholesteric liquid in 
the presence of an electric field is similar. The condi- 
tion dg/dc = fixes c since, for a given a, fc or c must 
satisfy: 



2£{k) 



for a > ac ~ 2/1: 



(12) 



The period of (j){u) modulation is T = 2fc lC{k) in units 
of ^. For each such period along the helical pitch axis, (j) 
increases by tt, the director unwinds once, and therefore 
one loses 1 of the T^qo/ir = 2ak /C(fc)/7r twists imprinted 
over the interval T. The imprinting efficiency is given by 
the fractional number of twists lost: 



Co 



2akK.{k)-TT 
2aklC{k) 



(13) 



Provided that eo is close to unity, much of the chirality 
can be preserved. 

For small chiral power a < ac = 2/7r, the minimising 
condition dg/dc = has no solution and the minimum 
free energy occurs exactly 'At c — 1. The cusp has the 
energy density g = (2/7r)^ — 1, and a logarithmically di- 
vergent period T which implies eo = 1. The director gets 
arrested, attempting to untwist but never actually man- 
ages a full turn within a finite distance. The imprinting 
is therefore successful. 

For chiral power a > ac — 2/7r, minima are found for 
c > 1. The period is no longer infinite and twists are 
lost. For large a, i.e. the nematic penetration depth 
large compared to the cholesteric pitch, we can expand 
the elliptic function for small k to find a~c(l — — 
3fc''/128...). Alternatively k ~ l/a(l - l/'ia^...) and we 
find a period of T — 2k K,{k) ^ vr/a in units of ^. Thus 
for every actual distance of tt^/q; = n / q^, (j) accumulates 
an increment of tt. That is, cq — > 0, corresponding even- 
tually to the case of complete unwinding of the imprinted 
helical pattern. 
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FIG. 3. Imprinting efficiency vs. network parameter a. The 
slope at the critical point a — q+ diverges and the efficiency 
decays rapidly. 
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In writing equation (|^), we have assumed that the chi- 
ral solvent was completely replaced with an achiral one 
as was the case in experiment |^]. However, we can triv- 
ially generalise F to the case where the current solvent 
is chiral: 

F = j dz \ [K2 {0' - qf + sin2(0 - q^z)] , (14) 

where q is the wave number that the solvent would in- 
troduce in the absence of crosslinks. Our previous case 
corresponds to q = Q. The mathematical procedure re- 
mains applicable, with an adjustment of the parameter 
a = ^go to: 

a = S.{qo - q) 

which indicates a chiral environment of the same handed- 
ness can help preserve the chirality by reducing the value 
of a (negative values of a are mathematically equivalent 
to —a, if we switch the sign of the (j) modulation). A 
chiral environment of the other handedness, q < Q, in- 
stead increases a and has the opposite effect. It reduces 
imprinting, as one might expect. By simply varying 5, 
via the solvent composition, we have a powerful way to 
map out the wind-unwinding transition. 

Another interesting case is the induction of a 
cholcsteric state in an initially uniform ncmatic network 
upon introducing a chiral solvent |l0| . That is go = 
and a' = ^q, where q is the chiral pitch wave number due 
to the chiral solvent. In this case, eg given in Fig. 2 can 
be interpreted as the network resistance to imprinting - 
until the point a' reaches 2/tt there will be no helicity 
induced in the network at all, and thereafter it rises with 
a'. 

Finally, by varying the amount of solvent in the crossed 
network, we can also tune, i.e. contract or expand, the 
volume relative to its original: V — > f3V. With this, we 
expect from molecular interpretations, K2 — > K2 for a 
special case where the local nematic order is preserved 
(the nematic order variation due to the solvent or tem- 
perature will be examined elsewhere), but Di Di/f3 
as it depends on crosslink density, and finally assuming 
an isotropic expansion/contraction, qo — » qoP^^^"^. Thus 



we have 

so swelling can, albeit weakly, increase the parameter a 
in our model and discourage the preservation of the im- 
printed chirality. 

In conclusion, we have proposed a continuum model 
for chiral imprinting in nematic elastomers. The model 
predicts the residual chirality when the chiral solvent is 
removed. A twist-untwist transition emerges from the 
theory, and an experimental verification of this transi- 
tion would be a useful test of the theory presented here. 
Further work is focussed on the mechanical properties 
of networks; we expect an external mechanical field will 
have substantial infiuence on the imprinting of chirality. 
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